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ABSTRACT 

We investigate the impact of neutral hydrogen (H i) in galaxies on the statistics of 
21-cni fluctuations using analytic and semi-numerical modelling. Following the reioni- 
sation of hydrogen the H i content of the Universe is dominated by damped absorption 
systems (DLAs) , with a cosmic density in H i that is observed to be constant at a level 
equal to ~ 2% of the cosmic baryon density from z ~ 1 to z ~ 5. We show that ex- 
trapolation of this constant fraction into the reionisation epoch results in a reduction 
of 10-20% in the amplitude of 21-cm fluctuations over a range of spatial scales. The 
assumption of a different percentage during the reionisation era results in a propor- 
tional change in the 21-cm fluctuation amplitude. We find that consideration of Hi in 
galaxies/DLAs reduces the prominence of the H ii region induced shoulder in the 21-cm 
power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. 
We also estimate the 21cm-galaxy cross PS, and show that the cross PS changes sign 
on scales corresponding to the H ii regions. From consideration of the sensitivity for 
forthcoming low- frequency arrays we find that the effects of HI in galaxies/DLAs on 
the statistics of 21-cm fluctuations will be significant with respect to the precision of a 
PS or cross PS measurement. In addition, since overdense regions are reionised first we 
demonstrate that the cross-correlation between galaxies and 21-cm emission changes 
sign at the end of the reionisation era, providing an alternative avenue to pinpoint 
the end of reionisation. The sum of our analysis indicates that the H i content of the 
galaxies that reionise the universe will need to be considered in detailed modelling of 
the 21-cm intensity PS in order to correctly interpret measurements from forthcoming 
low- frequency arrays. 

Key words: cosmology: diffuse radiation, large scale structure, theory - galaxies: 
High redshift, inter-galactic medium 



1 INTRODUCTION 

The process of hydrogen reionisation is thought to have 
started with ionised (Hll) regions around the first galax- 
ies, which later grew to surround groups of galaxies. Reion- 
isation completed once these Hll regions overlapped and 
occupied most of the volume between galaxies. Much re- 
cent theoretical attention has focused on the power spec- 
trum (PS) of 21-cm emission from neutral hydrogen (Hi) 
during the reionisation era (Zaldarriaga, Furlanetto & Hern- 
quist 2004; Furlanetto et al. 2004; Morales et al. 2005; Bow- 
man et al. 2006). In particular, the ionisation structure of 
the intergalactic medium (IGM) owing to UV emission as- 
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sociated with star formation has been studied in detail us- 
ing analytic (Furlanetto et al. 2004; Barkana 2007), numer- 
ical (McQuinn et al. 2006; Iliev et al. 2008), and more re- 
cently, semi-numerical models (Zahn et al. 2007; Mesinger 
& Furlanetto 2007; Geil & Wyithe 2008). 

These studies describe a scenario in which very large 
Hll regions form around clustered sources within overdense 
regions of the IGM. The formation of these Hll regions has 
a significant efi'ect on the shape of the 21-cm PS because 
information on the small scale features in the density field 
is erased from the signal originating within the ionised re- 
gions. Conversely, the creation of large ionised regions im- 
prints large scale features on the distribution of 21-cm in- 
tensity. The sum of these effects is to move power from small 
to large scales, leaving a shoulder shaped feature on the PS 
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at the characteristic scale of the H ii regions (Furlanetto et 
al. 2004). The detailed morphology of the ionisation struc- 
ture will therefore yield information about both the ionising 
sources and the structure of absorbers in the IGM on small 
scales (McQuinn et al. 2006). As reionisation leaves a strong 
imprint on the 21-cm PS, measuring the latter has become a 
key goal for learning about the reionisation epoch (Furlan- 
etto, Oh & Briggs 2006, and references therein). In addi- 
tion, since reionisation is driven by galaxy formation, whose 
statistics reflect those of the underlying density field, the 
detection of the redshifted 21-cm signal will not only probe 
the astrophysics of reionisation, but also the matter PS dur- 
ing the epoch of reionisation (McQuinn et al. 2006; Bow- 
man, Morales & Hewitt 2007). Finally, it has been recognised 
that cross-correlating galaxy surveys with 21-cm maps could 
yield powerful new insights into the morphology of reioni- 
sation, as well as eliminate some of the difficulties related 
to foreground removal (Furlanetto & Lidz 2007; Wyithe & 
Loeb 2007; Lidz et al 2008). With these ideas as motiva- 
tion, several experiments are currently under development 
that aim to detect the 21-cm signal during reionisation, in- 
cluding the Low Frequency Array^ (LOFAR) , the Murchison 
Widefield Array^ (MWA) and the Precision Array to Probe 
Epoch of Reionization^ (PAPER), and more ambitious de- 
signs are being planned such as the Square Kilometer Array* 
(SKA). 

Up until now only the component of H I residing in the 
IGM has been considered in relation to forecasts of the star 
tistical 21-cm signal. However, after the completion of reion- 
isation there is known to be a residual Hi fraction of a few 

percent in high density clumps which are believed to reside 
within galaxies (e.g., Prochaska et al 2005). This high den- 
sity contribution to the H I content of the Universe is also 
present during the reionisation era (we refer to this high 
density Hi as galactic Hi throughout the paper). Moreover, 
because galaxies at high redshift are biased relative to the 
density field, this galactic Hi contribution could provide a 
significant perturbation to the predicted statistics of 21-cm 
fluctuations. 

Wyithe & Loeb (2007) have modelled the density de- 
pendent reionisation process using a semi-analytic model 
that incorporates the important physical processes associ- 
ated with galaxy bias and radiative feedback. In agreement 
with rmmerical simulations (McQuinn et al. 2006; Ilicv et 
al. 2008), this model demonstrates that galaxy bias leads to 
enhanced reionisation in overdense regions. In this paper we 
use this model and its semi-numerical extension to explore 
the effect that the H I content of the biased galactic sources 
of reionisation would have on the fluctuations in redshifted 
21-cm emission. We begin by describing our density depen- 
dent model of reionisation in § 2. We then summarise the ob- 
served evolution of the H I density in § 3, and our results for 
the evolution of the 21-cm signal in § 4. We next describe the 
effect of galactic Hi on the 21-cm signal using analytic and 
semi-numerical models in § 5-7, and discuss prospects for 
detection in § 8. We conclude in § 9. Throughout this paper 
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we adopt a concordance cosmology for a flat ACDM uni- 
verse, (nm,nA,nb,fe,cr8,rt) = (0.27,0.73,0.046,0.7,0.8,1), 
consistent with the constraints from Komatsu et al. (2008). 
All distances are in comoving units unless stated otherwise. 



2 DENSITY DEPENDENT ANALYTIC 
MODEL OF REIONISATION 

In regions of the IGM that are overdense, galaxies will be 
over-abundant for two reasons: first because there is more 
material per unit volume to make galaxies, and second be- 
cause small scale fiuctuations need to be of lower amplitude 
to form a galaxy when embedded in a larger scale ovcrdcn- 
sity (the so-called galaxy bias; see Mo & White 1996). Re- 
garding reionisation of the IGM, the first effect will result 
in a larger density of ionising sources. However this larger 
density will be compensated by the increased density of gas 
to be ionised, which also increases the recombination rate. 
The process of reionisation also contains several layers of 
feedback. Radiative feedback heats the IGM and results in 
the suppression of low-mass galaxy formation (Efstathiou, 
1992; Thoul & Weinberg 1996; Quirm et al. 1996; Dijk- 
stra et al. 2004). Such feedback effects can potentially be 
more intense in overdense regions, leading to weaker galaxy 
formation bias than might be expected from simple linear 
bias models (Kramer et al 2006). Probing the morphology 
of reionisation could potentially lead to constraints on such 
feedback effects. 

To compute the relation between the local overdensity 
and the brightness temperature of redshifted 21-cm emission 
we use the model described in Wyithe & Loeb (2007). Here 
we summarise the main features of the model, and describe 
the additions made for the purpose of including the possible 
contribution from galaxies. The evolution of the ionisation 
fraction by mass Qs,r of a particular region of IGM with 
scale R and overdensity 5 (at observed redshift Zohs) may be 
written as 



dt 



0.76 



dF,oiiS,R,z,Mio 



+ (1 - QS.R. 



dt 

dFcol{S, R,Z,Mmln) 



dt 



V D{Zohs)J 



(1) 



where Mon is the number of photons entering the IGM per 
baryon in galaxies, as is the case-B recombination coeffi- 
cient, C = 2 is the clumping factor (which we assume, for 
simplicity, to be constant), and D{z) is the growth factor 
between redshift z and the present time. The production 
rate of ionising photons in neutral regions is assumed to be 
proportional to the collapsed fraction Fcoi of mass in ha- 
los above the minimum thresholds in neutral (Mmin), and 
in ionised (A/ion) regions. We assume Mmin to correspond 
to a virial temperature of 10*K, representing the hydrogen 
cooling threshold, and A/ion to correspond to a virial tem- 
perature of lO^K, representing the mass below which infall 
is suppressed from an ionised IGM (Dijkstra et al. 2004). 
In a region of co-moving radius R and mean overdensity 
S{z) = 5D{z)/D{zohs) [specified at redshift z instead of the 
usual a = 0], the relevant collapsed fraction is obtained 
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from the extended Press-Schechter (1974) model (Bond et 
al. 1991) as 



Fcoi{S,R,z) = evic\ 



-5{z) 



(2) 



where erfc(a;) is the complementary error function, a^{R) 
is the variance of the density field smoothed on a scale R, 
and (Tgai is the variance of the density field smoothed on a 
scale -Rgai, corresponding to a mass scale of Mmin or Mion 
(both evaluated at rcdshift z rather than at 2: = 0). In this 
expression, the critical linear overdensity for the collapse of 
a spherical top-hat density perturbation is 5c « 1.69. 

The model assumes that on largo (linear regime) scales 
most ionising photons are absorbed locally, so that the ion- 
isation of a region is caused by nearby ionisation sources. 
This assumption is certainly justified during the eaxly stages 
of reionisation, when the mean free path for ionising pho- 
tons is short. However, oven later in the reionisation process, 
the mean free path always remains smaller than the charac- 
teristic H II bubble size (it could be smaller if mini-halos or 
pockets of residual H I block ionising photons between the 
sources and the edge of the H 11 region.) Our local ionisation 
assumption is therefore valid as long as the characteristic 
bubble size is smaller than the spatial scale of the correla- 
tions we consider. This requirement is met in regimes where 
the fraction of regions at a particular scale are fully ionised 
is very low. This point is discussed quantitatively in Wyithe 
& Morales (2007). 



3 THE H I CONTENT OF GALAXIES 

Following the completion of reionisation, around 2% of the 
hydrogen content of the Universe is observed to be in Hi 
within galaxies. Moreover, this two percent fraction is ob- 
served to be roughly constant from 2; ~ 1 to z ~ 5 (e.g. 
Prochaska et al. 2005). However, to estimate the effect of 
galactic Hi on 21-cm fluctuation statistics during the reion- 
isation epoch we must extrapolate to higher redshift. There 
are two competing factors. On the one hand the collapsed 
fraction of mass rises as reionisation progresses, which would 
imply a galactic H I fraction that decreases towards high red- 
shift. On the other hand, the fact that the Hi fraction is 
constant below « ~ 5 suggests that the fraction of hydro- 
gen within galaxies that is in Hi increases with redshift. In 
this paper we conservatively assume that the fraction of gas 
within galaxies that is H I increases at a rate that preserves 
the 2% galactic H I mass fraction (relative to the total hy- 
drogen content of the Universe) into the reionisation epoch. 
At redshifts where there is not a sufficient collapsed fraction 
to maintain the 2% value, the galactic H I fraction is allowed 
to drop below 2%. 

This is obviously a fairly simplistic model. The actual 
physics of the HI content in galaxies depends both on star 
formation and feedback processes within the galaxy, as well 
as on the intergalactic UV background. For instance, as the 
UV background falls toward high redshift, the HI fraction 
within galaxies could increase, and perhaps conspires to can- 
cel with the decreasing collapse fraction to maintain the ob- 
served roughly constant Hi fraction of ~ 2%, as seen from 
2 ~ 1 — 5. Indeed, the fraction in protogalactic Hi may be 



even larger than the ~ 2% we assume in the pre-reionisation 
epoch, particularly if there is a significant population of 
minihalos^ with Tvir < lO'^K. Modelling this quantity is 
difficult and beyond the scope of this paper. For now, we 
note that our empirical extrapolation becomes increasingly 
uncertain at high redshift; it is probably reasonable in the 
redshift range 2; ~ 6 — 8, where some of the most interesting 
effects (such as the change of sign of the cross-correlation 
between galax;ies and 21-cm emission) take place. 



4 EVOLUTION OF THE 21-CM INTENSITY 
SIGNAL 

To illustrate the contribution of galactic H I, we use the fol- 
lowing example. We first find the value of iVion that yields 
overlap of ionised regions at the mean density IGM by ~ 6 
(Fan et al. 2006; Gnedin & Fan 2006; White et al. 2003), and 

then integrate equation (1) as a function of 5 and R. At a 
specified redshift, this yields the filling fraction of ionised 
regions within the IGM on various scales J? as a function of 
overdensity. We may then calculate the corresponding 21- 
cm brightness temperature contrast owing to Hi in both 
the IGM and galaxies 



T{5, R) = 22mK 



1 + z 
7.5 



1/2 



(1 - Q5,H - Q^r) (1 + ^-5 



(3) 



where the pre- factor of 4/3 on the overdensity refers to the 
spherically averaged enhancement of the brightness tem- 
perature due to peculiar velocities in overdense regions 
(Bharadwaj & Ah 2005; Barkana & Loeb 2005), and 1-Qs°k 
is the mass averaged fraction of H I in galaxies relative to the 
total hydrogen content of the Universe. Note that we only 
account for brightness temperature fluctuations due to vari- 
ation in the density and ionisation fraction of the IGM. We 
ignore variations in the spin temperature (which we assume 
to be always much greater than the CMB temperature) due 
to fluctuations in heating and radiative Lya coupling; these 
are important primarily early in the reionisation process, 
when galaxies are in any case rare. 



4.1 The auto-correlation function of 21-cm 
intensity fluctuations 

Since the underlying probability distribution of overdensi- 
ties dP/d6 is known (a Gaussian of variance aD), we may 
compute the observed probability distribution for T, 



dP dP 
dT^^^'^lS 



dS_ 
dT 



(4) 



The second moment of this distribution corresponds to 
the auto-correlation function of brightness temperature 



^ Since these are photoevaporated as reionisation progresses 
(Shapiro ct al 2004), they will reside primarily in neutral regions, 
and present an additional biased Hi contribution that is anti- 
correlated with galaxies but which we do not model here. 



© 0000 RAS, MNRAS 000, 000-000 



4 Wyithe, Warszawski, Geil, Oh 





z 



Figure 1. Modcl-I associating galactic Hi witli star forming galaxies. Upper Left: The evolution of the mean brightness temperature with 
redsliift. Upper Right: The evolution of the mean Hi fraction with redshift. The observed points are inferred from observations of damped 
hya absorbers (Prochaska et al. 2005). Lower Left: The evolution of the auto-correlation function of brightness temperature with redshift 
(computed within regions of size 10'). Lower Right: The evolution of the cross-correlation function with redshift (computed within regions 
of size 10'). In each case the grey and dark lines correspond to models which exclude and include the galactic HI respectively. The small 
upper sections of each panel show the difference between the models without galactic Hi and including galactic Hi. 



smoothed on an angular radius 6: 



on an angular radius 6 is 



1 

.V2na{R) 

where 

(T> 



l-^^j d5{T{6,R)-{T)fe-^ 



/ d5 T(5, R) 

2na{R) J ^ ' 



, (5) 



(6) 



and (t(-R) is the variance of the density field (at redshift z) 
smoothed on a scale R. 



1 

V^a{R) 



I 



dS {T{S, R) 



,(7) 



where Sg — bS, and b{Mg,z) is the galaxy bias for a halo 
of mass Mg at redshift z (Shcth, Mo & Tormcn 2001). To 
calculate the cross-correlation function we assume Mg = 
1O"M0 for the host mass of observed galaxies throughout 
this paper. This host mass yields a density of galaxies that 
is comparable to the density observed in the Subaru Deep 
Fields (Kashikawa et al. 2006) at « ~ 6.6 (Furlanetto & 
Lidz 2007). 



4.2 The cross-correlation between galaxies and 
21-cm intensity fluctuations 

The inclusion of galactic H I in the calculation of the 21-cm 
signal will also modify predictions for the cross-correlation of 
21-cm emission with galaxies. The cross-correlation function 
of galaxy overdensity and brightness temperature smoothed 



5 MODEL-I: GALACTIC H I ASSOCIATED 
WITH STAR FORMING GALAXIES 

In this section we use our analytic model to estimate the 
effect of galactic Hi on 21-cm fluctuations. We assume that 
the star forming galaxies responsible for reionisation are the 
harbourers of galactic H I, and calculate the fraction of galac- 
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Figure 2. Model-I associating galactic Hi with star forming galaxies. Left: The angular dependence of the auto-correlation function of 
brightness temperature (computed at z = 7). Right: The angular dependence of the cross-correlation function (computed at z = 7). In 
each case the grey and dark lines correspond to models which exclude and include the galactic H I respectively. 



tic hydrogen that is in Hi form 



/hi.coI = 
min ( 1, 



0.02 



F,^l{z, Mn,in) (1 - 0) + F,,l{z, Mion)Q 



(8) 



where Fcoi and Q are computed for the mean universe {5 — 0, 
R = oo). This enforces the condition that the fraction of HI 
in galaxies is either ~ 2%, or equal to the (lower) collapse 
fraction in galaxies. We can then compute the fraction of 
cosmic hydrogen that is galactic H I in regions of overdensity 
S and radius R 



Qcol 
6,R 



/hi,co1 [Qs,RFcol{d, R, Z, Mion) 

+ (1 - Q6,r) Fcoi{6, R, z, M^i„)] . (9) 



This contribution to the H I fraction of the Universe can be 

used with equation (14) to estimate the effect of galactic Hi 
on the statistics of 21-cm emission. 

The evolution of mean brightness temperature, together 
with the mean Hi fraction {xm} = 1 — Q — Q'^°^ are plot- 
ted in Figure 1. The data points show the observed fraction 
from damped Lya absorbers for comparison (Prochaska ct 
al 2005). The solid dark lines show the case including the 
galactic H I contribution. The thick grey line shows the IGM- 
only model for comparison. The residuals plotted above each 
panel show the magnitude of the effect when the galaxy con- 
tribution is ignored. The inclusion of galactic Hi increases 
the global 21-cm signal slightly as expected. 

Examples of auto-correlation functions are shown in 
Figures 1 and 2. The auto-correlation function is plotted 
as a function of redshift at 6 = 10' (Figure 1), and as a 
function ol 6 aX z = 7 (Figure 2). In Figure 1 the residuals 
are plotted above each panel to show the magnitude of the 
error introduced when the galaxy contribution is ignored. 
Early in the reionisation process, before the appearance of 
Hll regions begins to dominate the fluctuation amplitude, 
the inclusion of a galactic fraction enhances the 21-cm fluc- 
tuations. However, in contrast to the mean signal, the size of 
fluctuations are reduced by the presence of galaxies late in 
the reionisation era. This reduction can be traced to the fact 
that galaxies are biased towards overdense regions, while 
neutral IGM is biased towards underdense regions which 
are reionised last. As a result, the inclusion of a galactic 
H I fraction reduces the intensity contrast between overdense 
and underdense regions. The size of the residuals shows that 



within this model, the 2% galactic H I fraction reduces the 
fluctuation amplitude over a large fraction of the reionisa- 
tion epoch. At ^; ~ 7, it reduces the fluctuation amplitude by 
around 10% at all angles (which makes sense, since galactic 
HI constitutes about ~ 10% of the total HI content of the 
universe at that point). 

Examples of the cross-correlation functions axe shown 
in Figure 1 as a function of redshift at fixed angle {6 = 10'), 
and in Figure 2 as a function of angle at fixed redshift 
(z = 7). The inclusion of galactic Hi reduces the amplitude 
of the cross-correlation function by about 10% (relative to 
the case where galactic Hi is ignored) late in the reioni- 
sation era for the reasons discussed in § 4.2. In addition, 
the cross- correlation changes sign at overlap because the H I 
content of the Universe shifts from being dominated by the 
underdense regions of IGM to H I in galaxies (which reside in 
overdense regions). This sign change will provide an unam- 
biguous pointer to the redshift at which reionisation ends. 
Moreover, because the sign change will occur over a narrow 
frequency interval, it should provide an important check on 
possible sources of systematic error that could be present in 
measurements of the PS of intensity fluctuations which is 
always positive by construction. 



6 MODEL-II: VARYING HOST MASS FOR 
GALACTIC H i 

The calculations presented in Figures 1 and 2 assume that 
the galactic Hi contributing to the 21-cm emission is found 
in the star forming systems. However it is possible that the 
connection between galactic H I and star formation is not di- 
rect, either because there is Hi remaining in older low mass 
galaxies, or because the fraction of galactic gas that is Hi 
is host mass dependent. As presented, our model is not able 
to address either of these issues. Therefore, to ascertain the 
possible range of influence that galactic Hi might have on 
the 21-cm fluctuation statistics, we have assumed a charac- 
teristic mass Mcoi for the hosts of galactic H I, and replaced 
the expression for Tigm (equation 14) with 



TiGM {5, R) 



22mK 



l + z 
7.5 



1/2 



(1-Q«,«) ( 1+ ^-^j +a;Hi,coi (l+ - 



,(10) 
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Figure 3. Model-II with fixed characteristic host mass for galactic H I. Left: The evolution of the auto-correlation function of brightness 
temperature with redshift (computed within regions of size 10'). Right: The evolution of the cross-correlation function with redshift 
(computed within regions of size 10'). In each case the grey and dark lines correspond to models which exclude and include galactic Hi 
respectively. The small upper section of each panel shows the difference between the models without and including gala<;tic H I. 




where we have defined the mass averaged fraction of cosmic 
hydrogen that is galactic Hi (a::Hi,coi). We note that in or- 
der to maintain self consistency, the host mgiss of the galac- 
tic Hi must correspond to a collapsed fraction [Fcoi(Mcoi)] 
that is in excess of a::Hi,coi. At a particular redshift the max- 
imum possible host mass (Mcoi,x) corresponding to a col- 
lapsed fraction that equals xhi.coI is given by evaluation of 

Fcoi(5,i?, 2,Mcoi,x) =a;ni.coi. (11) 

The host mass at which the galaxy bias (6) in equation (10) 
is evaluated is therefore 

M = min (Mcoi.x, Mcoi) . (12) 

With these modifications for the galactic Hi distribution, 
we may compute the 21-cm fluctuation statistics as before 
using equation (7). 

The resulting auto-corrclation functions are shown in 
Figure 3 as a function of redshift at fixed angle {d = 10'). 
Three values for Mcoi are shown, IO^Mq, 1O^°M0 and 
10^^ Mq. Late in the reionisation era, the amplitude of the 
auto-correlation function is decreased by ~ 10% (relative to 
the case where galactic H I is ignored) . Below 2 ~ 6 larger. 



more biased hosts of galactic H I result in a more significant 
increase of the fluctuation amplitude. However at z > 6 
Mcoi,x < Mc, and hence our model predicts the same evolu- 
tion for different values of Mcoi (see equation 12). The cross- 
correlation function is also shown in Figure 3. The cross- 
correlation changes sign at overlap as before. More massive 
hosts of the galactic H I result in a cross-correlation with a 
larger amplitude following reionisation. 



7 SEMI-NUMERICAL SIMULATIONS 

The analytic calculations described thus far provide a quali- 
tative description of the effect of galactic H I on the statistics 
of 21-cm fluctuations and their association with galaxies. 
However, this analytic model is unable to describe fluctua- 
tions on scales comparable to the characteristic bubble scale 
once reionisation becomes established (Wyithe & Morales 
2007). In the next section of this paper we therefore inves- 
tigate the statistics of 21-cm fluctuations including galac- 
tic Hi in a semi- numerical model for the reionisation of a 
three-dimensional volume of the IGM (Geil & Wyithe 2008; 
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Zahn et al. 2007; Mesinger & Furlanetto 2007). Our semi- 
numerical simulations are based on model-II with respect to 
the galactic H I contribution. 

7.1 The ionisation field 

Our modelling follows the procedure outlined in Geil & 
Wyithe (2008), and we refer the reader to that paper for 
details of the model. The model employs a semi-analytic 
prescription for the reionisation process, which is combined 
with a realisation of the density field. We construct an ionisa- 
tion field based on a Gaussian random field for the overden- 
sity of mass, combined with the value of the ionised fraction 
Qs,R (equation 1) as a function of overdensity 5 and smooth- 
ing scale R. Wc repeatedly filter the linear density field at 
logarithmic intervals on scales comparable to the box size 
down to the grid scale size. For all filter scales, the ionisa- 
tion state of each grid position is determined using Qs^r and 
deemed to be fully ionised if Q6,r ^ 1. All voxels within a 
sphere of radius R centred on these positions are fiagged and 
assigned Qs,r ~ 1, while the remaining non-ionised voxels 
are assigned an ionised fraction of Qs^Rf^^^, where i?f,min 
corresponds to the smallest smoothing scale. A voxel forms 
part of an H II region if Qs.r > 1 on any scale R. In this 
paper we present simulations corresponding to a linear den- 
sity field of resolution 256^ , with a comoving side length of 
512 Mpc. This procedure yields an ionisation map Qion(x) 
as a function of position x. 

Having computed the ionisation field we then find the 
distribution of the galactic Hi component. Rather than try 
to assign individual galaxies to the simulation we make the 
approximation that Poisson noise will be negligible and as- 
sign a smooth density of H I 

Phi,coi(5) = a;Hi,coiPH[l -|- bd{x)], (13) 

where pni.coi and pn are the mass averaged densities of Hi 
that is collapsed in galaxies, and the total hydrogen density 
respectively. The position dependent brightness temperature 
of the simulation box becomes 

TiGM(f) = 22mKf^^j 

X [(1 - Qion) (1 + 3)+ xmMl + bS)] . (14) 

Note that our semi-numerical model does not compute pecu- 
liar velocities, and so equation (14) does not include a pecu- 
liar velocity induced enhancement of the brightness temper- 
ature in overdense regions. As in earlier sections we consider 
a model in which the mean IGM is roionised at z = 6. Wo 
again assume that star formation proceeds in halos above 
the hydrogen cooling threshold in neutral regions of IGM. 
In ionised regions of the IGM star formation is assumed to 
be suppressed by radiative feedback (see §2). 

7.2 Variance in 21-cm emission 

Since our semi-numerical model computes the three dimen- 
sional ionisation structure in the IGM, including the effect 
of H II regions, we can use it to compute the evolution of the 
variance. The resulting auto-correlation functions are shown 
in the left hand panel of Figure 4 as a function of redshift at 
fixed angle {6 = 10'). As in the analytic model presented in 



Figure 3, three values for Mcoi are shown, 10" Mo, IO^'Mq 
and 10" Mq. These curves can be compared directly with 
the analytic approximation. This comparison shows that 
while the analytic model yields the correct qualitative be- 
haviour, it is not quantitatively correct, both in terms of 
the redshift and amplitude of maximum fluctuations. How- 
ever, in agreement with the analytic calculation, our semi- 
numerical calculations of the auto-correlation function show 
that the fluctuation amplitude is modifled at the level of 
10-20% (relative to the case where galactic Hi is ignored) 
by the presence of a 2% galactic H I fraction throughout 
the reionisation era. The semi-numerical calculations of the 
cross-correlation function are also shown in the right hand 
panel of Figure 4 as a function of redshift at fixed angle 
[9 = 10'). Comparison with Figure 3 also shows that the 
analytic calculation of the cross-correlation is qualitatively 
correct, but does not predict the correct quantitative evolu- 
tion. 



7.3 21-cm power spectrum 

Figure 5 shows 21-cm PS computed from the semi- numerical 
simulations at z = 7 (left) and 2 = 9 (right). We plot the 
dimensionless PS A2i(fc) = fc^/(27r^)P2i(fc), where P21 is 
the PS of 21-cm fluctuations. Both the case of the IGM 
alone (A|i), and the case assuming a host mass of Mcoi = 
10^° Mq for the galactic Hi (A21 coi) shown (lower pan- 
els). We also show the relative fluctuation (52i,coi = (A^i — 
Aii.coO/Aii) of the 21-cm PS owing to galactic Hi (upper 
panels) . 

These figures illustrate the shoulder in the PS that cor- 
responds to the typical bubble scale and which is due to the 
movement of power from small to large scales that accom- 
panies the formation of Hii regions. Because the galactic 
H I is biased towards H 11 regions, this movement of power 
is lessened when the contribution of galactic Hi is consid- 
ered. Hence the PS at largo scales evaluated from simulations 
which include H I in gala^cies is lower than the PS computed 
using the IGM alone. The effect of galactic H I is therefore 
to change the shape of the 21-cm PS rather than just the 
amplitude. 



7.4 21cm-gal£ixy cross power spectrum 

The central and upper panels of Figure 6 show the modulus 

of the dimensionless 21cm-galaxy cross PS with (A2i_g_(.oi) 
and without (A^i^g) galactic H I as well as the galactic H I in- 
duced fiuctuation [(52i,g,coi = (Aii^g — Aji^g^coO/Aii.g], com- 
puted from the semi-numerical simulations at 2; = 7 (left) 
and z = 9 (right). Again a host mass of Mcoi = W^°Mq is 
considered for the galactic Hi. An observed galaxy mass of 
Mg = 10^'^ Mq is assumed. In the lower panels of Figure 6 
we show the corresponding coefficient of the 21cm-galaxy 
cross PS. 

Since overdense regions, where galaxies are concen- 
trated, are reionised first these figures show an anti- 
correlation on large scales, which drops in strength to zero 
on small scales after the formation of H II regions. We find 
that the inclusion of galactic H I lessens the amplitude of the 
anti-correlation, since a fraction of H I is now co-located with 
the galaxies inside the H 11 regions. In addition, the cross PS 
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Figure 5. Lower panels: 21-cm power spectra computed using scmi-numcrical models at z = 7 (left) and 2 = 9 (right). In each case 
models are shown for a host mass of M^^i = IO^^Mq for the galactic Hi, as well as for the IGM-only case. Upper panels: The relative 
perturbation on the PS due to galactic Hi [<52i,col = (^ii ~^2i col)/^2i]- Also plotted for comparison is the estimated sensitivity (lower 
panels), and the sensitivity relative to A|j^ (upper panels) within bins of width Afe = k/10 assuming 1000 hr integration on 1 field for 
the MWA and MWA5000 (dotted and thin-grey lines). 



changes sign on small scales which reflects the correlation of 
galaxies with the galactic H I inside the H ii regions, where 
no power is contributed in 21-cm fluctuations of the IGM. 
Thus, with the caveat that galaxies can be selected without 

bias from IGM absorption (which may not be true if, for in- 
stance, the galaxies are selected in Lya emission — see Lidz 
et al 2008 for discussion for the latter case), the scale at 
which the 21cm-galaxy cross PS changes sign could be used 
to probe the scale of H ii regions late in the reionisation era. 



8 SENSITIVITY TO THE EFFECT OF 

GALACTIC Hi ON 21-CM FLUCTUATIONS 

Before concluding, we compute the sensitivity with which 
the effect of galactic H I could be detected using forthcoming 
low-frequency arrays. 



8.1 Sensitivity to the 21-cm PS 

To compute the sensitivity AP21 (k) of a radio- interferometer 
to the 21-cm PS, we follow the procedure outlined by Mc- 
Quinn et al. (2006) and Bowman, Morales & Hewitt (2007) 
[see also Wyithe, Loeb & Geil (2008)]. The important issues 
are discussed below, but the reader is referred to these pa- 
pers for further details. The sensitivity to the PS comprises 
components due to the thermal noise, and due to sample 
variance within the finite volume of the observations. We 
consider observational parameters corresponding to the de- 
sign specifications of the MWA, and of a hypothetical fol- 
lowup to the MWA (termed the MWA5000). In particular 
the MWA is assumed to comprise a phased array of 500 
tiles. Each tile contains 16 cross-dipoles to yield an effec- 
tive collecting area of Ac — 16(A^/4) (the area is capped for 
A > 2.1m). The physical area of a tile is Auic = 16m^. The 
tiles are distributed according to a radial antenna density 
of p(r) oc r~^, within a diameter of 1.5km and outside of a 
flat density core of radius 18m. The MWA5000 is assumed 
to follow the basic design of the MWA. The quantitative 



differences are that the telescope is assumed to have 5000 
tiles within a diameter of 2km, with a flat density core of 
80m. In each case we assume 1 field is observed for 1000 
hr. Following the work of McQuinn et al. (2006) we as- 
sume that foregrounds can be removed over 8MHz bins, 
within a bandpass of 32MHz [foreground removal there- 
fore imposes a minimum on the wave-number accessible of 
fcn,in ~ 0.04[(1 + 2)/7.5]~iMpc-i]. 

The sensitivity to the 21-cm PS per mode may be writ- 
ten 



sP2i{k,e) 



D^AD / 
Ai^tint n{k±) \Ae 



sky 



+ P2i{k,e), (15) 



where D is the co-moving distance to the centre of the survey 
volume which has a co-moving depth AD. Here n{k±) is 
the density of baselines which observe a wave vector with 
transverse component k±, and 6 is the angle between the 
mode k and the line of sight. The thermal noise component 
(first term) is proportional to the sky temperature, where 
Tsky ~ 250 (iy^)^ *^ K at the frequencies of interest. The 
second term corresponds to sample variance. The overall 
sensitivity is 



AP2i{k,e) ^ 5P2i{k,e)/y^N^{k,e), 



(16) 



where Nc{k,9) denotes the number of modes observed in a 
fe-space volume d^k (only modes whose line-of-sight com- 
ponents fit within the observed band-pass axe included). 
In terms of the fc-vector components k and 6, Nc = 
2TTk'^smedkdeV/{2TTf where V = D'^AD{\^/Ae) is the ob- 
served volume. Taking the spherical average over bins of 9, 
the sensitivity to the 21-cm PS is 

' E r.o (17) 



[AP2i(fc)]2 



[AP2i(k,e)]^ 



The spherically averaged sensitivity curves for the 
MWA (within bins of Ak = fc/10) are plotted as the dotted 
lines in the lower panels of Figure 5. The sensitivity as a 
ratio of the PS computed without a galactic Hi contribu- 
tion {AP21/P21) is plotted in the upper panels of Figure 5 
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Figure 6. 21cm-galaxy cross PS computed using semi-numerical models. Central and Lower panels: The modulus and coefRcient of the 

21cm-galaxy cross PS respectively computed at z = 7 (left) and 2 = 9 (right). In each case models are shown assuming a host mass of 
-'^col = 10""' Mq for the galactic Hi, as well as the IGM-only case. The host mass of observed galaxies is assumed to be Mg = IO'^^Mq. 
?7pper panels: The relative perturbation on the PS due to galactic Hi [<52l,g,col = (■^21 g ~^2l g col)/'^2l gl- ^1^° plotted for comparison 
is the estimated sensitivity (central panels) , and the sensitivity relative to A| ^ (upper panels) within bins of width Afc = A:/ 10 assuming 
1000 hr integration on 1 field for the MWA and MWA5000 (dotted and thin-grey lines). To estimate this sensitivity galaxies were assumed 
to be observed down to Mg = 10^^ Mq over the entire field. 



(again as dotted lines). These estimates illustrate that late 
in reionisation, the effect of the galactic H I will be at a level 
above the sensitivity of the first generation low-frequency ar- 
rays. However, the low resolution of the MWA would mean 
that the shape-change of the PS owing to the galactic Hi 
would not be detected until very close to overlap. The cor- 
responding sensitivity curves for the MWA5000 are plotted 
as the thin grey lines in Figure 5. The larger collecting area 
of the MWA5000 would allow the effect of galactic H I to be 
detected out to higher redshifts and at smaller scales. 

8.2 Sensitivity to the 21cm-galaxy cross PS 

To compute the sensitivity of the MWA to the 21cm-galax;y 
cross PS (P2i,g) we follow the discussion of Purlanetto & 
Lidz (2007). The sensitivity to a particular mode is 

2[5P2i,^{k,e)f = \P2i.Jk,e)f +5P2i{k,e)5P^{k,9), (18) 

where <5Pg(A;, 9) = b^Ps{k, 9) +11^^ is the uncertainty in the 
galaxy PS, and Pe is the underlying mass PS® at redshift 

® Note that for consistency with our estimate of the 21-cm PS, 
we are neglecting redshift space distortions in calculation of the 
galaxy PS. 



z. In the second term rig is the density of galaxies (which 
we approximate as rig = Mgdn/dM^, where dn/dM is the 
Press-Schechter 1976) mass function and Mg is the halo mass 
of galaxies in the survey^. After calculation of the total sen- 
sitivity 

AP2i,g(A;, 9) = SP2i,s{k, 9)/y/Nc{k,9), (19) 

and taking the spherical average over bins of 9 as before, the 
sensitivity to the 21cm-galaxy cross PS is 

[AP^iAkW =^ [AP2iAk,eW 

The spherically averaged sensitivity curves for the 
MWA (within bins of Afc = fc/10) are plotted as the dotted 
lines in the central panels of Figure 6. The sensitivity as a 
ratio of the fiducial IGM only PS (AP2i,g/f'2i,g) is also plot- 
ted in the upper panels of Figure 6. Here we have assumed 
that galaxies have been detected down to a resolution limit 
corresponding to a host mass of Mg = 10^^ M© over the full 
MWA field. These sensitivity estimates illustrate that the 
effect of the galactic Hi on the cross-correlation between 

^ Note that we have neglected redshift errors when computing 
the Poisson component of the uncertainty in the galaxy PS. 
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galaxies and 21-cm emission could be detected with the first 
generation low-frequency arrays. Moreover the resolution of 
the MWA would be sufficient to detect the sign-change of 
the 21cm-galax;y cross PS at the scale of the Hll regions 
late in reionisation. The corresponding sensitivity curves for 
the MWA5000 are plotted as the thin grey lines in Figure 6. 
The larger collecting area of the MWA5000 would greatly 
increase the significance with which the effect of galactic H I 
could be measured. 



9 CONCLUSION 

In this paper we have investigated the impact of Hi in 
galaxies on the statistics of 21-cm fluctuations using ana- 
lytic and semi-numerical models. Our models are unable to 
self-consistently compute the galactic H I content of galax;ies 
prior to the end of reionisation. As an input to our model we 
have therefore assumed that during the reionisation era 2% 
of hydrogen is in the form of H I and located within galax- 
ies. This number is motivated by observations of the mass 
weighted fraction of cosmic hydrogen in Hi after the end 
of reionisation, which is dominated by damped absorption 
systems and which has a constant value of ~ 2% between 
2 ~ 1 and z ^ 5. Our modelling shows that this assump- 
tion results in a reduction of 10-20% in the amplitude of 
21-cm fluctuations over a range of spatial scales, and over 
a large fraction of the reionisation era. In addition to the 
amplitude of 21-cm fluctuations we have also modelled the 
cross-correlation between galaxies and 21-cm emission. We 
find that the inclusion of galaxies decreases the amplitude of 
the cross-correlation by a few tens of percent. In addition, 
the cross-correlation between galaxies and 21-cm emission 
will change sign at the end of the reionisation era, providing 
an alternative avenue to pinpoint the end of reionisation. 

We find that our analytic estimates become less applica- 
ble once H ii regions become a dominant feature of the IGM. 
We have therefore supplemented our analytic estimates with 
semi- numerical modelling of the three dimensional ionisation 
structure of the IGM. We find that the qualitative conclu- 
sions from our analytic calculations are vindicated by this 
modelling. In addition to the above estimates of auto and 
cross-correlation, we have used our semi-numerical model to 
compute 21-cm power spectra, and 21cm-galaxy cross power 
spectra. The effect of galactic Hi is to change the shape of 
the 21-cm PS rather than just the amplitude. In particular, 
because the galactic Hi is biased towards Hll regions, we 
find that the H II region induced shoulder in the PS is less 
prominent when the contribution of galactic H I is consid- 
ered. 

The amplitude of the 21-cm PS is modified in a scale de- 
pendent way by up to 20% when a 2% galactic H I fraction is 
included. Of course the galactic H I fraction is very uncertain 
at high redshift. However we find that the fractional modi- 
fication of the 21-cm fluctuation statistics is approximately 
proportional to the density of galactic Hi. Therefore, if we 
assume a value lower(higher) than the 2% observed at z < 5, 
then our results for the error introduced through neglect 
of galactic Hi will be over(under)estimated. On the other 
hand, the change of sign in the cross-correlation between 
galaxies and the 21-cm signal is robust to our assumptions. 

We find that the inclusion of galactic Hi lessens the 



amplitude of the anti-correlation between galaxies and 21- 
cm emission since a fraction of Hi is now co-located with 
the galaxies inside the Hll regions. In addition, the cross PS 
changes sign on small scales, which refiects the correlation of 
galaxies with the galactic H I inside the H ii regions, where 
no power is contributed in 21-cm fluctuations of the IGM. 
Thus, the scale at which the 21cm-galaxy cross PS changes 
sign could be used to probe the scale of H II regions late in 
the reionisation era. 

We have estimated the sensitivity of the MWA to the 
spherically averaged 21-cm PS and 21cm-galaxy cross PS. 
Our calculations illustrate that the effect of the galactic H I 
on 21-cm fluctuations is at a level that would be significant 
with respect to the sensitivity of the MWA late in reionisa- 
tion. However the low resolution of the MWA would mean 
that the shape change of the PS owing to the galactic H I 
could not be detected until very close to overlap. On the 
other hand, when combined with a suitable galaxy redshift 
survey, the resolution of the MWA would be sufficient to de- 
tect the sign change of the 21cm-galaxy cross PS at the scale 
of the H II regions late in reionisation. A foUowup telescope 
comprising 10 times the collecting area of the MWA would 
measure the effects of galactic H I with high significance. 

In summary, our modelling shows that the H I content 
of the galaxies that reionise the universe provides a signifi- 
cant contribution to the statistics of 21-cm fluctuations. The 
galactic H I contribution to the 21-cm intensity will therefore 
need to be considered in detailed modelling of the 21-cm in- 
tensity PS in order to correctly interpret measurements from 
the next generation of low-frequency arrays. 
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